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Abstract 

We discuss the baseline optical configuration for the Laser Interferometer Space Antenna (LISA) 
mission, in which the lasers are not free-running, but rather one of them is used as the main fre- 
quency reference generator (the master) and the remaining five as slaves, these being phase-locked 
to the master (the master-slave configuration). Under the condition that the frequency fluctuations 
due to the optical transponders can be made negligible with respect to the secondary LISA noise 
sources (mainly proof- mass and shot noises), we show that the entire space of interferometric com- 
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also be constructed with the master-slave system design. The corresponding hardware trade-off 
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needed for implementing the one-way method, and the master-slave configuration, are essentially 
O ■ identical. Either operational mode could therefore be implemented without major implications on 

^ ' the hardware configuration. 

bX). 

We then derive the required accuracies of armlength knowledge, time synchronization of the 

^ ■ onboard clocks, sampling times and time-shifts needed for effectively implementing Time-Delay 

■ 

Interferometry for LISA. We find that an armlength accuracy of about 16 meters, a synchronization 
accuracy of about 50 ns, and the time jitter due to a presently existing space qualified clock will 
allow the suppression of the frequency fluctuations of the lasers below to the level identified by 
the secondary noise sources. A new procedure for sampling the data in such a way to avoid the 
problem of having time shifts that are not integer multiples of the sampling time is also introduced, 
addressing one of the concerns about the implementation of Time-Delay Interferometry. 
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I. INTRODUCTION 



LISA, the Laser Interferometer Space Antenna, is a three-spacecraft deep space mis- 
sion, jointly proposed to the National Aeronautics and Space Administration (NASA) and 
the European Space Agency (ESA). The LISA scientific objective is to detect and study 
low-frequency cosmic gravitational radiation by observing phase differences of laser beams 
interchanged between drag- free spacecraft J]. 

Modeling each spacecraft with two optical benches, carrying independent lasers, frequency 
generators (called Ultra Stable Oscillators), beam splitters and photo receivers, the measured 
eighteen time series of frequency shifts (six obtained from the six one-way laser beams 
between spacecraft pairs, six from the beams between the two optical benches on each of 
the three spacecraft, and six more from modulation of the laser beams with USO data) 
were previously analyzed by Tinto et al. [2j. There it was shown that there exist several 
combinations of these eighteen observables which exactly cancel the otherwise overwhelming 
phase noise of the lasers, the phase fluctuations due to the non-inertial motions of the six 
optical benches, and the phase fluctuations introduced by the three Ultra Stable Oscillators 
into the heterodyned measurements, while leaving effects due to passing gravitational waves. 

The analysis presented in j^J relied on the assumptions that (i) the frequency offsets of any 
pair of independent lasers (assumed there to be ~ 300 MHz) could be observed within the 
detection bandwidths of the photo receivers where the one-way Doppler measurements are 
performed, and (ii) the telemetry data rate needed by two of the three spacecraft to transmit 
their measured one-way Doppler data to the third spacecraft (where the interferometric 
combinations are synthesized) is adequate. Although the technology LISA will be able 
to use should make possible the implementation of Time-Delay Interferometry (TDI) as 
discussed in |2|, the possibility of optimizing the design of the optical layout, while at the 
same time minimizing the number of Doppler data needed for constructing the entire space 
of interferometric observables, was not analyzed there. Here we extend those results to a 
different optical configuration, in which one of the six lasers is the provider of the frequency 
reference (albeit time-delayed) for the other five via phase-locking. This master-slave optical 
design could provide potential advantages, such as smaller frequency offsets between beams 
from pairs of different lasers, hardware redundancy, reliability, and can result in a smaller 
number of measured data. An outline of the paper is given below. 
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In Section II we summarize TDI, the data processing technique needed for removing 
the frequency fluctuations of the six lasers used by LISA, and other noises. In order to 
show that the entire space of interferometric observables LISA can generate can also be 
reconstructed by using a master-slave optical configuration, we consider the simple case of 
spacecraft that are stationary with respect to each other. After showing that the entire 
space of interferometric observables can be obtained by properly combining four generators, 
(a, /3,7, £), we then derive the expressions for the four generators corresponding to the 
master-slave optical configuration. By imposing some of the one-way measurements entering 
into (a, P, 7, C) to be zero (the so called locking conditions) , we show that the expressions for 
these generators can be written in terms of the one-way and two-way Doppler measurements 
corresponding to the locking configuration we analyzed. Section III and Appendix A provide 
a theoretical derivation and estimation of the magnitude of the phase noise expected to be 
generated by an optical transponder. In Section IV we analyze and compare the hardware 
requirements needed for implementing both optical designs, while in Section V we turn to 
the estimation of the armlength and time synchronization accuracies, as well as time-shift 
and sampling time precisions needed for successfully implementing TDI with LISA. Our 
comments and conclusions are finally presented in Section VI. 

II. TIME-DELAY INTERFEROMETRY 

The description of TDI for LISA is greatly simplified if we adopt the notation shown 
in Figure 1, where the overall geometry of the LISA detector is defined. The spacecraft 
are labeled 1, 2, 3 and distances between pairs of spacecraft are L±, L 2 , L 3 , with Li being 
opposite spacecraft i. Unit vectors between spacecraft are n«, oriented as indicated in figure 
1. We similarly index the phase difference data to be analyzed: s 3 \ is the phase difference 
time series measured at reception at spacecraft 1 with transmission from spacecraft 2 (along 
L 3 ). Similarly, s 2 i is the phase difference series derived from reception at spacecraft 1 with 
transmission from spacecraft 3. The other four one-way phase difference time series from 
signals exchanged between the spacecraft are obtained by cyclic permutation of the indices: 
1 — > 2 — > 3 — ^ 1. We also adopt a useful notation for delayed data streams: 531,2 = ssi(t— L2), 
531,23 = ss\(t — L 2 — L3) = 531,32, etc. (we take the speed of light c = 1 for the analysis). 
Six more phase difference series result from laser beams exchanged between adjacent optical 
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FIG. 1: Schematic LISA configuration. Each spacecraft is equidistant from the point O, in the plane 
of the spacecraft. Unit vectors hi point between spacecraft pairs with the indicated orientation. 
At each vertex spacecraft there are two optical benches (denoted 1, 1*, etc.), as indicated. 

benches within each spacecraft; these are similarly indexed as nj = 1, 2, 3 ; i ^ j). 

The proof-mass-plus-optical-bench assemblies for LISA spacecraft number 1 are shown 
schematically in figure 2. We take the left-hand optical bench to be bench number 1, while 
the right-hand bench is 1*. The photo receivers that generate the data S21, S31, T21, and T31 
at spacecraft 1 are shown. The phase fluctuations of the laser on optical bench 1 is pi(t); 
on optical bench 1* it is p\{t) and these are independent (the lasers are for the moment 
not "locked" to each other, and both are referenced to their own independent frequency 
stabilizing device). We extend the cyclic terminology so that at vertex % (i = 1,2,3) the 
random displacement vectors of the two proof masses are respectively denoted Si(t) and S*(t), 
and the random displacements (perhaps several orders of magnitude greater) of their optical 
benches are correspondingly denoted Aj(t) and A*(t). As pointed out in Q], the analysis does 
not assume that pairs of optical benches are rigidly connected, i.e. Aj 7^ A*, in general. The 
present LISA design shows optical fibers transmitting signals both ways between adjacent 
benches. We ignore time-delay effects for these signals and will simply denote by fii(t) the 
phase fluctuations upon transmission through the fibers of the laser beams with frequencies 
Ui, and v*. The n%(t) phase shifts within a given spacecraft might not be the same for 
large frequency differences V{ — v*. For the envisioned frequency differences (a few hundred 
megahertz), however, the remaining fluctuations due to the optical fiber can be neglected j^J. 
It is also assumed that the phase noise added by the fibers is independent of the direction 
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of light propagation through them. 

Figure 2 endeavors to make the detailed light paths for these observations clear. An 
outgoing light beam transmitted to a distant spacecraft is routed from the laser on the local 
optical bench using mirrors and beam splitters; this beam does not interact with the local 
proof mass. Conversely, an incoming light beam from a distant spacecraft is bounced off the 
local proof mass before being reflected onto the photo receiver where it is mixed with light 
from the laser on that same optical bench. The inter-spacecraft phase data are denoted S31 
and s 2 i m figure 2. Beams between adjacent optical benches within a single spacecraft are 
bounced off proof masses in the opposite way. Light to be transmitted from the laser on an 
optical bench is first bounced off the proof mass it encloses and then directed to the other 
optical bench. Upon reception it does not interact with the proof mass there, but is directly 
mixed with local laser light, and again down converted. These data are denoted T31 and T21 
in figure 2. 

The terms in the following equations for the s»j and phase measurements can now 
be developed from figures 1 and 2, and they are for the particular LISA configuration in 
which all the lasers have the same nominal frequency Uq, and the spacecraft are stationary 
with respect to each other. The analysis covering the configuration with lasers of different 
frequencies and spacecraft moving relative to each other was done in Q], and we refer the 
reader to that paper. 

Consider the S3i(t) process (equation 0) below. The photo receiver on the left bench of 
spacecraft 1, which (in the spacecraft frame) experiences a time- varying displacement Ai, 
measures the phase difference S31 by first mixing the beam from the distant optical bench 
2* in direction n 3 , and laser phase noise p\ and optical bench motion Ag that have been 
delayed by propagation along L3, after one bounce off the proof mass (Si), with the local 
laser light (with phase noise pi). Since for this simplified configuration no frequency offsets 
are present, there is of course no need for any heterodyne conversion. 

In equation (j3J) the T31 measurement results from light originating at the right-bench laser 
(p\, AJ), bounced once off the right proof mass [o~D, and directed through the fiber (incurring 
phase shift /ii(t)), to the left bench, where it is mixed with laser light (pi). Similarly the 
right bench records the phase differences S21 and Ti\. The laser noises, the gravitational 
wave signals, the optical path noises, and proof-mass and bench noises, enter into the four 
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FIG. 2: Schematic diagram of proof-masses-plus-optical-benches for a LISA spacecraft. The left- 
hand bench, 1, reads out the phase signals S31 and T31. The right hand bench, 1*, analogously 
reads out S21 and T21. The random displacements of the two proof masses and two optical benches 
are indicated (lower case 5i for the proof masses, upper case Aj for the optical benches.) 

data streams recorded at vertex 1 according to the following expressions 
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Eight other relations, for the readouts at vertices 2 and 3, are given by cyclic permutation 
of the indices in equations (QJ-(|1J). 

The gravitational wave phase signal components, sf™ , — 1, 2, 3, in equations (HJ and 
(JHJ) are given by integrating with respect to time the equations (1), (2) of reference |5[ that 
relate metric perturbations to frequency shifts. The optical path phase noise contributions, 
^opt. path^ w j 1 j c j 1 includes shot noise from the low signal-to-noise ratio (SNR) in the links 
between the distant spacecraft, can be derived from the corresponding term given in 0. 
The Tij measurements will be made with high SNR so that for them the shot noise is 
negligible. 

The laser-noise-free combinations of phase data can readily be obtained from those given 
in [2I for frequency data. We use the same notations: X, Y, Z, a, (3, 7, (, etc., but the 
reader should keep in mind that here these are phase measurements. 

The phase fluctuations, s^- , 7y , i,j = 1, 2, 3, are the fundamental measurements needed 
to synthesize all the interferometric observables unaffected by laser and optical bench noises. 
If we assume for the moment these phase measurements to be continuous functions of time, 
the three armlengths to be perfectly known and constant, and the three clocks onboard 
the spacecraft to be perfectly synchronized, then it is possible to cancel out exactly the 
phase fluctuations due to the six lasers and six optical benches by properly time-shifting 
and linearly combining the twelve measurements s^- , Ty , i,j = 1,2,3. The simplest 
such combination, the totally symmetrized Sagnac response (, uses all the data of Figure 2 
symmetrically 

C = s 32,2 — -323,3 + s 13,3 — s 31,l + s 21,l — $12,2 

+ 2i( r 23 ~~ r 13),12 + ( r 31 — r 2l),23 + ( r 12 ~ r 32),13] 

+ ^[(T23 - r 13 ),3 + (^31 - T 2 l),l + ("712 - T 32 ),2), (5) 

and its transfer functions to instrumental noises and gravitational waves are given in [^] 
a nd B I respectively. In particular, { has a » 8i x-pu. 8 e response" to gravitationa. radiation, 
i.e. a 5-function gravitational wave signal produces six distinct pulses in £ |5J, which are 
located with relative times depending on the arrival direction of the wave and the detector 
configuration. 

Together with (, three more interferometric combinations, (a,(3,j), jointly generate the 
entire space of interferometric combinations 0, Q]. Their expressions in terms of the 
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measurements s,-,-, follows 

Oi = S21 — S31 + S 132 — Si2,3 + 332,12 — S 2 3 5 13 + - [(t 23 — T 13 ) )2 + (t" 23 — T 13 ) 13 + (r 31 — T 21 ) 
+ (r 3 i - r 2 i),i23 + (T12 - T"32),3 + (r 12 - r 32 ), 12] , (6) 

with and 7 derived by permuting the spacecraft indices in a. Like in the case of (, a 
^-function gravitational wave produces six pulses in a, /3, and 7. In equations (0 EJ) it is 
important to notice that the Ty measurements from each spacecraft always enter into the 
interferometric measurements as differences taken at the same time. This property naturally 
suggests a locking configuration that makes these differences equal to zero, as we will show 
in the next section. 

We remind the reader that the four interferometric responses (a, (3, 7, () satisfy the fol- 
lowing relationship 

C - C.123 = a ,l - « 23 + /3,2 - Al3 + 7,3 ~ 7,12 • (7) 

Jointly they also give the expressions of the interferometric combinations derived in jjj, j^: 
the Unequal-arm Michelson (X, Y, Z), the Beacon (P, Q, R), the Monitor (E, F, G), and the 
Relay (U, V, W) responses 

X i = a 23 - P,2 ~ 7,3 + C , (8) 
P = C-«i, (9) 
E = a - C,i , (10) 
U = 7,i-0> (11) 

with the remaining expressions obtained from equations (jHl El E3 HH) by permutation of the 
spacecraft indices. All these interferometric combinations have been shown to add robustness 
to the mission with respect to failures of subsystems, and potential design, implementation, 
or cost advantages 0,0. 



A. Locking Conditions 

The space of all possible interferometric combinations can be generated by properly time 
shifting and linearly combining the four combinations (a, j3, 7, (), as given above. Although 
they have been derived by applying TDI to the twelve one-way Doppler data, in what follows 
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FIG. 3: Simplified optical layout of the LISA interferometer, showing all the optical benches, proof 
masses and lasers. 

we will show that they can also be written in terms of properly selected and time shifted 
two-way and one-way Doppler measurements. These can be generated by phase locking five 
of the six lasers to one of them, as it is described below. 

Assume, without loss of generality, the laser on bench 1* to be the master. Although 
there are several other possible locking schemes, the one chosen minimizes the number of 
locking conditions between the master and any given slave. Furthermore, locking schemes 
relying on more than one master could be implemented, but we will not address those in 
this paper. We will also assume the spacecraft to be stationary relative to each other. This 
assumption simplifies the analysis, and does not affect the validity of the general result sj]. 
Under this assumption, the frequency provided by the master laser 1* can be used as input 
reference for the slaved lasers. In other words the slaves will then have the same center 
frequency as the master, and their phase fluctuations will be related to the fluctuations of 
the master laser as well as any other fluctuations introduced into the received light beam 
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prior to reception and locking. 

In order to understand the topology of the beams as the various slaves are locked to the 
master, let us follow the light paths from the master laser, 1*, to the slaves, as shown in 
Figure 3. Let us start first with the light beam that is bounced off the back of the proof-mass 
1*. This beam is then directed to the other bench, 1, where it is used as the input frequency 
reference for the laser there, and the measurement T31 is made. Light is then re-transmitted 
back to bench 1* where the measurement T21 is performed. Since the phase of the laser 1 is 
locked to that of the master, the relative phase fluctuations r 31 can be adjusted as follows 

731 = 7-21 , (12) 

where we have assumed the noise introduced by phase-locking to be negligible (on this point 
see the theoretical derivations in Section III and Appendix A). Similarly, light beams from 
lasers 1 and 1* are transmitted to the lasers on the benches 2* and 3 respectively, where 
the lasers are locked to the incoming beams. From benches 2* and 3 beams are transmitted 
to the lasers on benches 2 and 3* respectively, where again locking is performed similarly 
to what is done onboard spacecraft 1. Finally, along arm 1, only two one-way relative 
phase measurements can be performed as it is easy to see. Since for the moment we have 
assumed a configuration with stationary spacecraft, the optical configuration described above 
can be translated into the following locking conditions on some relative phase fluctuation 
measurements 

7-31=7-21 , 7-13 = 7-23 , 7"32 = T12 , S 2 3 = S 32 = . (13) 

The locking conditions define specific relationships among the phase fluctuations from vari- 
ous noise sources. As an example, the conditions 

7-31 = t 2 i , s 2 3 = , (14) 

imply the following relationships among the laser phase fluctuations, the proof-mass noises, 
the gravitational wave signal, and the two bench noises onboard spacecraft 1 

Pi = Pi + "a n 3 • {81 - Ai) + u n 2 ■ (<5* - A*) , (15) 
= s 23 + s 23*' Path + Pi, 2 -P3 + V0 ~ 2n 2 • 5 3 + n 2 ■ A 3 + n 2 ■ A* 2 , (16) 

with similar expressions following from the other locking conditions given in equation (|13p. 

11 



If we now substitute the locking conditions into the expressions for (a, ft, 7, (), we obtain 
their expressions in terms of the remaining measurements 

Clo. — s 13,3 ~~ -531,1 + -521,1 ~ -512,2 > (17) 

Olio. = S 2 l - S31 + -513,2 - Sl2,3 , (18) 

ftlo. = Sl2 + $21,3 + -513,23 ~ -531,12 , (19) 

llo. = S 13 - S 3 i, 2 + S 2 1,13 - -512,23 , (20) 

where the data (si 2 , s% 3 ) are one-way, and (s 2 i, s 3 i) are effectively two-way Doppler mea- 
surements due to locking. 

The verification that the combinations (ati ., ftlo., lio., Clo.) exactly cancel the laser phase 
fluctuations as well as the fluctuations due to the mechanical vibrations of the optical 
benches, can be performed by substituting the locked phase processes (such as [T3J fTHjl 
into the expressions for (si 2 , si 3 , s 2 i, s 3 i) (which are given by equations (JTJ EJ) and their 
permutations), and by further replacing them into equations (1171 - 12011 . 

The main result of implementing locking is quantitatively shown by equations (|17I-I2()|1 in 
that the number of measurements needed for constructing the entire space of interferometric 
combinations LISA will be able to generate is smaller by a factor of three than the number 
of measurements needed when only one-way data are used. 

Once (atio., ftio., lio., Cio.) are constructed according to the expressions given in equations 
()17I-I20j). all the other interferometric combinations can be derived by applying the identities 
given in equations (|H1- ITT]) . As an example, it is straightforward to show that equation (jHJ) 
implies the following expression for the unequal-arm Michelson combination Xi . 

Xio. = [S 2 1 - S3l] - [-521,33 - 531,22] , (21) 

which coincides with the expression for X, derived for the first time in j^J, in terms of the 
two-way Doppler measurements from the two LISA arms. 

As a final comment, we have analyzed also several locking configurations needed when 
the spacecraft are moving relative to each other. We have found that there exist techniques, 
when locking is implemented, which are similar to the one analyzed in for removing the 
noise of the onboard Ultra Stable Oscillators from the phase measurements. The conclusions 
derived above for the case of stationary spacecraft are therefore general, and an analysis 
covering locking configurations with moving spacecraft is available in 
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III. PHASE LOCKING PERFORMANCE 



The locking conditions given in equation (|13|) reflect the assumption that the noise due 
to the optical transponders is negligible. In this section we analyze the noise added by the 
process of locking the phase of the local laser to the phase of the received light. This noise 
will be in addition to the optical path and USO noises, which we consider separately. 

A block diagram of the phase locking control system is shown in figure BJa). The system 
consists of a photoreceiver, phasemeter, controller, actuator and laser (see section IV for 
a description of these subsystems). Each of these subsystems can be characterized by its 
transfer function and, when it applies, by a noise contribution jjj. The main inputs to the 
system are the phase noise of the local laser, Pl(s), and the phase fluctuations of the signal 
beam from the distant spacecraft, ps{s) (with s = a + iu being the Laplace variable). The 
closed loop output is the phase noise of the retransmitted laser beam, Pcl(s). 

From the block diagram shown in figure 4a, it is easy to see that the closed loop output 
phase Pcl{s) can be written in terms of the free-running laser phase noise, Pl(s), the input 
signal phase fluctuations, Ps{s), and the various feedback components' transfer functions 
and noises shown in Figure 4, as follows 

Pcl =p L + L{N A + AG [N M + M{N R + R {N D +p s - Pcl})}} ■ (22) 

This equation can be solved for p C L 

p L + L{N A + AG [N M + M(N R + R {N D + p s })}} 
PCL ~ 1 + LAGMR ' (23) 

where we have denoted with N R (s), N M (s), and N A (s) the noises due to the photoreceiver, 
phasemeter, and actuator respectively. The detection noise, Nd(s) is the error in the mea- 
surement of the relative phase of the two beams. This noise is the fundamental limit to the 
phase measurement process for the LISA detection system (see Appendix A). Since the link 
from the phasemeter to the controller will be digital, the controller noise will be due only 
to the finite precision of the digital phase information. We thus assume this noise to be 
negligible and do not include it in this analysis. 

Since (i) the product M(s) R(s) = 1 (the phase at the phasemeter output ought to be 
equal to the phase at the input point of the photoreceiver), (ii) the contribution of the 
actuator noise to the output, L(s)Na{s), is much smaller than the noise from the free- 
running laser (the free running laser noise is measured with the actuators attached and 
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FIG. 4: a) Block diagram of phase locking control system showing transfer functions and noise 
contributions of various system components, (b) Simplified control system block diagram. 

thus intrinsically contains this noise source), and (iii) the laser transfer function, L(s), is a 
passive low pass filter with a pole at several GHz (and so can be ignored in this discussion), 
we conclude that the block diagram in 0| (a) now simplifies to that shown in figure EJ^b), 
with the closed loop output now given by 

An [AT J_ 1 PL + AG [N T +ps] ( . 

Pcl = Pl + AG [N t + p s - Pcl] — ► Pcl = \ + AG ^ ' 

p CL - Ps ~N T + ^ (forAG»l). (25) 

The quantity Nt{s) is the total noise at the input to the controller and is given by 

N T = N M + MN R + N D ) . (26) 

Equation (|25|) shows that phase locking will drive the phase of the local laser, pcl(s), to 
that of the received laser, ps(s), with an error introduced by two terms. The first term is 
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the total measurement noise in the phase measurement, Nt(s), which is in turn determined 
by the detection noise, photoreceiver noise and phasemeter noise (see equation I2T?)). Recall 
that the detection noise, Njj(s), is the fundamental error in the measurement of the relative 
phase of the two beams for the LISA detection system. A rigorous calculation of this 
noise source, included in appendix |XJ shows that its root power spectral density is equal 
to ~ 1 /icycle/VSz at the output of the phase meter. Nr(s) is the electronic noise of the 
photoreceiver at the beat note frequency. Referenced to the output of the phasemeter this 
will be well below the 1 /icycle/vTIz level. A phasemeter noise floor of ~ 1 /xcycle/\/Hz has 
been set as a requirement for the phasemeter noise. This level of performance has already 
been demonstrated (^(J, over the frequency range of interest (1 mHz to 1 Hz) albeit 
with heterodyne frequencies of a few kilohertz . Given these estimates of the individual 
noise sources a total measurement error, Nj>(s) of less than 2 /^cycles/ vTfz is expected. 

The second source of error in Eq. (|25|) represents the finite suppression of the free running 
laser noise. This term is inversely proportional to the loop gain and so can be reduced by 
increasing the gain. Very high gains should be possible with LISA as the frequencies of 
interest are very low jjj, ~ 10 Hz and lower. A free- running laser frequency noise of 
1 MHz/-\/Hz at 1 mHz corresponds to a phase noise of 10 9 cycles/vTIz. A total loop gain 
of 10 15 is therefore required to suppress the contribution of laser frequency noise down to 
1 /icycle / \/Hz . At low frequencies the laser phase is altered by changing the temperature 
of the laser crystal. This actuator has a typical (voltage to frequency) gain of 5 GHz/Volt 
or 5 x 10 12 cycles/v^Hz at 1 mHz. Thus a controller gain of 200 V/cycle is needed. This 
gain requirement could be eased by "pre-stabilizing" the laser frequency, for example by 
locking to a low finesse cavity or other frequency reference. Initial results from bench top 
experiments (Ql,^^) indicate that loop gains of the order of 10 15 should be achievable, 
and hence that pre-stabilization may not be necessary. 

The master-slave configuration phase locking requirement should be that the noise intro- 
duced by the locking process is insignificant compared to the 20 /xcycle/ vTIz of optical path 
noise allocated in the LISA noise budget Q . As the analysis above has shown, phase locking 
is limited only by the measurement noise (assuming adequate gain). This measurement noise 
is common to both the one-way and the master- slave schemes and so, given adequate gain 
of the phase locking loop, there will be no difference in the performance of the two systems. 
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IV. HARDWARE REQUIREMENTS 



In this section we compare the hardware needed by the two implementations of TDI 
discussed in the previous sections. The discussion will focus on the minimum hardware 
requirements and will not fully consider redundancy or fallback options. We will consider 
LISA as composed of several basic subsystems, and compare the type and quantity of the 
components required by each scheme. 



A. Laser Frequency Stabilization System 

Both schemes rely on a sufficiently high laser frequency stability. This is because the 
cancellation of laser frequency noise via TDI is not exact due to finite accuracy of the arm 
length knowledge, finite timing accuracy, imperfect clock synchronization, and sampling 
time jitters. The frequency stabilization system could be composed of either an optical 
cavity, a gas cell, or a combination of both. The output is an error signal proportional 
to the difference between the laser frequency and the resonance frequency of the reference. 
Assuming that an optical cavity is used as the frequency reference with a Pound- Dr ever- Hall 
locking |l4j readout, the frequency stabilization system will consist of an electro-optic phase 
modulator, an optical cavity, a photoreceiver, a double-balanced mixer, and a low pass filter 
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B. Phasemeter 



A phasemeter is a device capable of measuring the phase of a photoreceiver output relative 
to the local USO. One could distinguish between two types of phasemeter: single-quadrature 
phasemeters, denoted PM(S), and full-range phasemeters, denoted PM(F). An example of a 
single-quadrature phasemeter is a mixer. It will have only a limited linear range as the output 
is generally sinusoidal. A single-quadrature phasemeter has several potential advantages over 
a full-range phasemeter including lower noise, higher speed operation, increased reliability 
and lower power consumption. However, it is restricted in usefulness to closed loop operation 
only and is therefore not suitable for use in the one-way method. A single-quadrature 
phasemeter could potentially be used in the master-slaves configuration where one laser is 
phase locked to another with a fixed phase shift (see section IV. E). 
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A full-range phasemeter has an output that is linearly proportional to phase over the 
entire range — n to ir. An example of a full-range phasemeter is a zero-crossing time interval 
analyzer. As we show in section V. E, phasemeters used in an open-loop configuration must 
have a dynamic range of at least ~ 10 10 at 0.1 mHz. 

C. Controller 

A controller takes the signal from either a frequency stabilization system or phasemeter, 
amplifies and filters it appropriately, and feeds it back to the frequency /phase actuators of 
a laser. A controller is needed to frequency lock a laser to a frequency reference or to phase 
lock one laser to another. In practice the frequency locking and phase locking controllers 
will differ by a pole in the controller transfer function and by a gain factor. Depending on 
the sophistication of its design, a controller could potentially be reconfigured in-flight to 
perform either function. 

D. Photoreceivers 

Each scheme requires twelve photoreceivers to measure the interference from the front and 
back of the proof masses. The photoreceiver unit will consist of a photodiode and low-noise 
electronic amplifiers. Although the photoreceivers for LISA will contain quadrant photodi- 
odes for alignment sensing, this is irrelevant for the following discussion which assumes that 
single element photodiodes are used. 

E. Requirements for the Master-Slave Configuration 

Table I summarizes the system components needed on each optical bench. The quantities 
shown represent the minimum requirements with no redundancy included. In the master- 
slave configuration scheme one master laser is frequency stabilized to its own frequency 
reference. All other lasers are phase locked in a chain to this master laser, as described 
in section IIA. For this reason the minimum system requirement is only one frequency 
stabilization system. However, the capability of stabilizing other lasers to a local stabilization 
system should be included for redundancy against failure of the master's stabilizing device. 
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Spacecraft /Bench 


Master-Slave Configuration 


One- Way Method 


SC 1* 


1 FS, 2 PR, 2 PM(F), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 


SC 1 


2 PR, 2 PM(F), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 


SC 2* 


2 PR, 1 PM(F), 1 PM(S), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 


SC 2 


2 PR, 2 PM(F), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 


SC 3* 


2 PR, 2 PM(F), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 


SC 3 


2 PR, 1 PM(F), 1 PM(S), 1 Controller 


1 FS, 2 PR, 2 PM(F), 1 Controller 



TABLE I: Comparison of minimum system components required for each scheme. FS: frequency 
stabilization system, PM(S): single-quadrature phasemeter, PM(F) full-range phasemeter, PR: 
photoreceiver. 



Providing each laser with a frequency stabilization system would provide a high level of 
redundancy and maintain compatibility with the one-way mode of operation. 

The master-slave configuration will require at least four full-range phasemeters for the 
main signal read out photoreceivers, s 2 i> s 3 i, Si 2 and Si 3 . Full-range phasemeters will also be 
required for measuring the signals derived from the back side of the proof masses, r^ , as the 
lasers on adjacent benches are locked by suppressing the difference of the phasemeter outputs, 
T ij~~ T kj — 0. If just one of these phasemeter outputs were used for phase locking then a single- 
quadrature phasemeter could be utilized. However, the noise in the optical fiber linking the 
benches, /ij, would then be imposed on the phase of the slave laser. Although this noise would 
be removed by including the second detector in the time-delay interferometry processing, 
this would impose unnecessary requirements on the stability of the fiber link to prevent 
increasing the slave laser's frequency noise. Furthermore, by suppressing the difference 
of these phasemeter outputs, we do not need to record this information for processing, 
as only the difference of the phasemeter outputs appears in the TDI equations (see for 
example, equations El and EJ). Single-quadrature phasemeters could potentially be used on 
the remaining two photoreceivers where phase locking of the beams returning to spacecraft 1 
(■S23 — S32 — 0) is performed. This is a relatively minor simplification, as suitable full-range 
phasemeters must be developed for the remaining ten photoreceivers. If the phase locking 
hierarchy must be reordered, for example due to a frequency stabilization system failure, 
then full-range phasemeters will also be required at these positions. Finally, implementing 
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Subsystem 


Master-Slave Configuration 


One- Way Method 


FS 


1 


6 


PR 


12 


12 


PM 


12 


12 


Controllers 


6 


6 


Observables 


4 


9 



TABLE II: Summary of minimum system components required for each scheme. FS: frequency 
stabilization system, PM: phasemeter, PR: photoreceiver, Observables: number of data streams 
required for processing. 

full-range phasemeters at all photoreceiver outputs will maintain compatibility with the one- 
way method. For these reasons we will drop the (F) or (S) suffix for the phasemeters and 
assume that only full-range phasemeters are used. 

Table II summarizes the total number of components required and the number of data 
streams to be recorded by each scheme. Using the master-slave configuration only four 
data streams remain to be measured («2i, S31, S12 and S13) as all other variables have been 
suppressed to effectively zero and therefore do not need to be recorded (although they should 
be monitored to ensure proper operation of the phase locking systems). One potential 
concern with the master-slave configuration is in the non-local nature of the control system, 
that is to say the main phase input to the control loop comes from the light from the 
distant spacecraft. The amplitude and phase of this beam could be adversely affected by 
many factors such as spacecraft alignment. Although this will also affect the quality of 
the one-way measurements, it could be more detrimental to the master-slave configuration. 
For example, the signal intensity could become so low to cause loss of phase lock entirely. 
Furthermore, because all the slaves are linked to the master by the phase-locking chain, if 
one phase- locking link is disrupted then all downstream links may also be lost. The severity 
of this non-local control problem depends on how often lock will be disrupted, and the 
difficulty of lock reacquisition. 
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F. Requirements for the One-way Method 

The one-way method employs a very symmetric configuration consisting of three identical 
spacecraft each containing two identical optical benches. The components needed on each 
optical bench are shown in the right hand column of Table I. The six lasers are frequency 
stabilized to their six respective frequency references. The phases of the beat notes of each 
local laser with the lasers from the adjacent spacecraft and bench are measured by full-range 
phasemeters to provide twelve data streams. The data from the phasemeters at the back of 
the proof masses on adjacent benches can be combined before being recorded without loss 
of generality. This reduces the total number of data streams to be recorded and processed 
to nine, as shown in Table II. However, as we will show in section V C the number of 
data streams to be exchanged between spacecraft will be the same for the one-way and 
master-slaves configuration. 

From Tables I and II it is clear that the one-way and master-slaves configurations are 
almost identical in terms of the quantity of components required. However, there are several 
more subtle differences in the hardware requirements of the two schemes. A disadvantage 
of the one-way method is that the laser frequencies may differ by as much as ±300 MHz 
if high finesse cavities are used as the frequency references. This maximum frequency off- 
set is determined by the free spectral range of the reference cavity, where a cavity with 
a round trip optical path of 0.5 m has been assumed. This large frequency offset will 
place greater demands on the photoreceivers' bandwidth than the master-slave configura- 
tion where frequency offsets can be kept to the minimum dictated by the Doppler shifts 
(less than ±10 MHz). Not only does the high heterodyne frequency place strict require- 
ments on the photodetector bandwidth, but also on the bandwidth stability. For example, 
assume that photoreceivers with a bandwidth of fi, w ~ 1 GHz are used for the main signal 
readouts. Although a heterodyne frequency, fh, of 300 MHz is within the 1 GHz photore- 
ceiver bandwidth there will be a 0.05 cycle phase delay at this frequency (the phase shift 
is equal to -aictan(fh/ fbw) radians for a simple single-pole type frequency response). If the 
bandwidth of a photodetector changes, then this phase shift will also change in a way that 
is indistinguishable from the effects of a gravitational wave. A simple calculation shows 
that a bandwidth change of a mere 0.023% (or 230kHz) would introduce a phase signal of 
10 /zcycles/ V / Hz for information at 300 MHz. If the heterodyne frequency is kept to 10 MHz 
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or less, then a photoreceiver bandwidth change of more than 0.6% (or 6 MHz) would be 
needed to produce a 10 /xcycles/vTIz phase shift. 

In section III it was shown that a total loop gain of 10 15 at 1 mHz is required to ensure 
that the phase locking loop performs correctly. Using the one-way method the phase locking 
loops are replaced by frequency locking to the reference cavity. The frequency stabilization 
system is expected to be limited by fluctuations in length of the cavity at a level of the 
order of 10 Hz/v^Hz at 1 mHz. Under this condition no advantage is gained by suppressing 
the measured laser noise below this level and so a loop gain of ~ 10 6 should suffice. The 
one-way method therefore requires a controller gain of only 2 x 10 -7 Volts/cycle at 1 mHz 
compared to the 200 Volts/cycle needed for the master-slave configuration. 



V. ACCURACY AND PRECISION REQUIREMENTS 

The limitations on the effectiveness of the TDI technique, either when the one-way or 
the master-slave configuration is implemented, come not only from all the secondary noise 
sources affecting the measurements Sjj , , i,j — 1,2,3 (such as proof-mass and optical 
path noises) but most importantly from the finite accuracy and precision of the quantities 
needed to synthesize the laser-noise free observables themselves. In order to synthesize the 
four generators of the space of all interferometric combinations, we need: 

(i) to know the distances between the three pairs of spacecraft; 

(ii) to synchronize the clocks onboard the three spacecraft, which are used in the data 
acquisition and digitization process; 

(iii) to be able to apply time-delays that are not integer multiples of the sampling time of 
the digitized phase measurements; 

(iv) to minimize the effects of the jitter of the sampling times themselves; and 

(v) to have sufficiently high dynamic range in the digitized data in order to be able to recover 
the gravitational wave signal after removing the laser noise. 

In the following subsections we will assume the secondary random processes to be due to 



the proof masses and the optical path noises 



We will estimate the minimum values of the 



accuracies and precisions of the physical quantities listed above that allow the suppression 
of the laser frequency fluctuations below the level identified by the secondary noise sources. 
For each physical quantity the estimate of the accuracy and/or precision needed will be 
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performed by assuming all the remaining errors to be equal to zero. Our estimates, therefore, 
will provide only an order of magnitude estimate of the accuracies and/or precisions needed 
for successfully implementing TDI. 

A. Armlength accuracy 

The TDI combinations described in the previous sections rely on the assumption of know- 
ing the armlengths sufficiently accurately to suppress laser noise well below other noises. 
Since the three armlengths will be known only within the accuracies 5 Li , i = 1,2,3 respec- 
tively, the cancellation of the laser frequency fluctuations from the combinations (a, j3, 7, £) 
will no longer be exact. In order to estimate the magnitude of the laser fluctuations remain- 
ing in these data sets, let us define Lj , i — 1, 2, 3 to be the estimated armlengths of LISA. 
They are related to the true armlengths Lj , i — 1, 2, 3, and the accuracies SLi , i — 1, 2, 3 
through the following expressions 

L l = L l + 5L l , i=l,2,3. (27) 

In what follows we will treat the three armlengths Lj , i — 1, 2, 3 as constants equal to 16.7 
light seconds. We will derive later on the time scale during which such an assumption is 
valid. We will also assume to know with infinite accuracies and precisions all the remaining 
physical quantities needed to successfully synthesize the TDI generators. 

If we now substitute equation (|2T|) into equations (0), and expand it to first order in 5 Li, 
it is easy to derive the following approximate expression for C(£), which now will show a 
non-zero contribution from the laser noises 

C(t) ~ ((t) + [P2,i3 - ^,12] SL X + [p 3>12 - p* 23 ] 6L 2 + [p 1>23 - p* 13 ] 6L 3 , (28) 

where the "'" denotes time derivative. Time-Delay Interferometry can be considered ef- 
fective if the magnitude of the remaining fluctuations from the lasers are smaller than the 
fluctuations due to the other noise sources entering in C(^); namely proof mass and optical 
path noises. This requirement implies a limit in the accuracies of the measured armlengths. 

Let us assume the six laser phase fluctuations to be uncorrelated to each other, their one- 
sided power spectral densities to be equal, the three armlengths to differ by a few percent, 
and the three armlength accuracies also to be equal. By requiring the magnitude of the 
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remaining laser noises to be smaller than the secondary noise sources, it is straightforward 
to derive, from Eq. ()28|) and the expressions for the proof mass and optical path noises 
entering into £(t) given in Pf, the following constraint on the common armlength accuracy 
\5L ( \ 



1 U sin 2 (7r/L) S p p ro ° f "(f) + ST lcal pat \f) 

l6Lcl ~ 277 V W) ■ (29) 

Here S p , S^ roof mas \ g^cai path are ^ e one _ s i(jed power spectral densities of the phase fluc- 
tuations of a stabilized laser, a single proof mass, and a single-link optical path respectively 
0| . If we take them to be equal to the following functions of the Fourier frequency / 0, [J 

S p (f) = 2.3 x 1CT 1 f- 8/s + 1.4 x 1(T 9 f- 27/5 cycles 2 Hz" 1 , (30) 

gproof mass^ = x 1Q -21 j-4 ? ^ 

goptical path^ = ^ x jq-10 cycleg 2 ^ ^ 

(where / is in Hz), we find that the right- hand- side of the inequality given by equation (|2U|) 
reaches its minimum of about 16 meters at the Fourier frequency f m i n = 1.0 x 10 -4 Hz, over 
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the assumed (10~ 4 , 1) Hz LISA band. This implies that, if the armlength knowledge \SL, 
can be made smaller than 16 meters, the magnitude of the residual laser noise affecting the 
( combination will be below that identified by the secondary noises. This reflects the fact 
that the armlength accuracy is a decreasing function of the frequency. For instance, at 10~ 3 
Hz the armlength accuracy goes up by almost an order of magnitude to about 155 meters. 

A perturbative analysis similar to the one described above can be performed for the 
remaining generators (a,/3, 7). We find that the correspondin&inequality for the armlength 
accuracy required for the a combination, \8L a \, is equal to |3j,|2| 

1 /[8 AnH^fL) + 16 si,iV/L)] SJ™" '"""(/) + 6 sf™ •"""(/) ,^ 

|,IJS Mf if*tf) ' (33) 

with similar inequalities also holding for (3 and 7. Equation (|33|) implies a minimum of the 
function on the right-hand-side equal to about 31 meters at the Fourier frequency f m i n = 
1.0 x 10~ 4 Hz, while at 10 -3 Hz the armlength accuracy goes up to 180 meters. 

Armlength accuracies significantly smaller than the level derived above can be achieved 
by implementing laser ranging measurements along the three LISA arms [ifl, and we do not 
expect this to be a limitation for TDI. 
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In relation to the accuracies derived above, it is interesting to calculate the time scales 
during which the armlengths will change by an amount equal to the accuracies themselves. 
This identifies the minimum time required before updating the armlength values in the TDI 
combinations. 

It has been calculated by Folkner et al. |l6| that the relative longitudinal speeds between 
the three pairs of spacecraft, during approximately the first year of the LISA mission, can 
be written in the following approximate form 



V h3 (t) = V$ sin (^j (ij) = (1,2) ; (1,3) ; (2,3) , (34) 



where we have denoted with (1,2), (1,3), (2,3) the three possible spacecraft pairs, vM* is a 
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constant velocity, and Tjj is the period for the pair In reference [3] it has also been 

shown that the LISA trajectory can be selected in such a way that two of the three arms' 
rates of change are essentially equal during the first year of the mission. Following reference 
Q, we will assume V$ = V$ ± V$ , with V$ = 1 m/s, V$ = 13 m/s, T lj2 = T lj3 w 4 
months, and T 2 3 rs 1 year. From equation ()34j) it is easy to derive the variation of each 
armlength, for example AL 3 (t), as a function of the time t and the time scale St during 
which it takes place 

r ( ) . / 2wt 



AL 3 (t) = V$> sin J St . (35) 

Equation ()35j) implies that a variation in armlength AL 3 w 10 m can take place during 
different time scales, depending on when during the mission this change takes place. For 
instance, if t <C T 12 we find that the armlength L 3 changes by more than its accuracy (~ 10 
meters) after a time St = 2.3 x 10 3 seconds. If however t ~ 7\ 2 /4, the armlength will change 
by the same amount after only St ~ 10 seconds instead. 



B. Clock synchronization accuracy 

The effectiveness of the TDI data combinations requires the clocks onboard the three 
spacecraft to be synchronized. Since the clocks will be synchronized with a finite accuracy, 
the laser noises will no longer cancel out exactly and the fraction of the laser frequency fluc- 
tuations that will remain into the TDI combinations will be proportional to the magnitude of 
the synchronization accuracy. In order to identify the minimum level of off-synchronization 
among the clocks that can be tolerated, we will proceed by treating one of the three clocks 
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(say the clock onboard spacecraft 1) as the master clock defining the time for LISA, and the 
other two to be synchronized to it. The relativistic (Sagnac) time-delay effect due to the 
fact that the LISA trajectory is a combination of two rotations, each with a period of one 
year, will have to be accounted for in the synchronization procedure. This is a procedure 
well known in the field of time-transfer, and we refer the reader to the appropriate literature 
for discussions on this point Here we will disregard this relativistic effect, and assume 
it can be compensated for with an accuracy better than the actual synchronization accuracy 
we derive below. 

Let us denote by 5t 2 , St^, the time accuracies (time-offsets) for the clocks onboard space- 
craft 2 and 3 respectively. If t is the time onboard spacecraft 1, then what is believed to be 
time t onboard spacecraft 2 and 3 is actually equal to the following times 

t 2 = t + St 2 , (36) 
h = t + 6t 3 . (37) 

If we now substitute equations (jHBllBTjl into the equation (j5J) for £, for instance, and expand 
it to first order in St^ , i = 2, 3, it is easy to derive the following approximate expression for 
C(t), which shows the following non-zero contribution from the laser noises 

C(t) ~ C(*) + [Pl,23 - Pl,l2 + P2.13 - P2.is] 8h + [p 2) 13 - #,23 + P3,12 ~ P3,ia] St 3 • (38) 

By requiring again the magnitude of the remaining fluctuations from the lasers to be smaller 
than the fluctuations due to the other (secondary) noise sources affecting ((t), it is possible 
to derive an upper limit for the accuracies of the synchronization of the clocks. If we assume 
again the six laser phase fluctuations to be uncorrelated to each other, their one-sided power 
spectral densities to be equal, the three armlengths to differ by a few percent, and the two 
time-offsets' magnitudes to be equal, by requiring the magnitude of the remaining laser noises 
to be smaller than the secondary noise sources it is easy to derive the following constraint 
on the time synchronization accuracy \8t^\ 
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with S p , SP roof mass , S° ptical path again as given in equations 

We find that the right-hand-side of the inequality given by equation (|3~5|l reaches its 
minimum of about 47 nanoseconds at the Fourier frequency f m in = 1-0 x 10~ 4 Hz. In other 
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words, clocks synchronized at a level of accuracy better than 47 nanoseconds will imply 
a residual laser noise that is smaller than the secondary noise sources entering into the ( 
combination. 

An analysis similar to the one described above can be performed for the remaining gen- 
erators (a, /3, 7). For them we find that the corresponding inequality for the accuracy in the 
synchronization of the clocks is now equal to 

J_ / [4 sin 2 (37r/L) + 8 sin^vr/L)] gg mas \f) + 3 jgCQ ^ 

' ^ " 27T/ V 4 S p(f) ' 

with equal expressions holding also for (3 and 7. The function on the right- hand- side of 
equation (J4T?|) has a minimum equal to 88 nanoseconds at the Fourier frequency f min = 
1.0 x 10 -4 Hz. As for the armlength accuracies, also the timing accuracy requirements 
become less stringent at higher frequencies. At 10 -3 Hz, for instance, the timing accuracy 
for ( and a, (3, 7 go up to 446 and 500 ns respectively. 

A 50 ns accuracy translates into a 15 meter armlength accuracy, which we argued earlier 
to be easily achievable by the use of laser ranging. We therefore expect the synchronization 
of the three clocks to be achievable at the level derived above. 



C. Telemetered signals and their sampling 

To reduce the LISA-to-Earth telemetry requirements, it is expected that the normal op- 
erational mode will not telemeter the phase time series to the ground directly. Rather, we 
expect the TDI observables to be computed at the LISA array, and then only the (relatively 
low data rate) laser-noise-free combinations transmitted to Earth. Thus, both implementa- 
tions of TDI discussed in this paper (the one-way method and the master-slave configuration) 
require phase measurements data to be exchanged among the spacecraft in order to synthe- 
size the four generators of the space of all interferometric combinations. Although it is clear 
that the master-slave configuration implies a smaller number of measurements than that 
required by the one-way method, the actual number of data that will need to be exchanged 
among the spacecraft can be made to be exactly the same for both, making their inter- 
spacecraft telemetry requirements identical. This can easily be understood by rewriting the 
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four generators (a, (3, 7, () in the following forms 

C = [sai.i - -531,1] + ^[( r 3i - r 2i),i + (V31 - r 2 i) j23 ] 

+ [-332,2 — s 12,2] + 2^ Tl2 ~~ T 32),2 + ( r 12 — r 32),13] 

+ [Sl3,3 - -523,3] + ^[(7"23 ~ T^s + (^23 ~ T 13 ) tl2 ] , (41) 

a = [s 21 - s 3 i] + ^ [(r 31 - t 2 \) + (r 31 - r 2 i),i 2 3] 

+ [^32,12 - Si 2i3 ] + - [(ri2 - r 32 ),3 + (Tl 2 - T 32 ),l 2 ] 

+ [-§13,2 - 823,13] + g t( T 23 _ T 13),2 + (T 2 3 - T13), 13] . (42) 

Equations (|4~T| l4*2*|) show that each generator can be formed by summing three different linear 
combinations of the data, each involving phase measurements performed onboard only a 
specific spacecraft. As an example, let us assume without loss of generality that ( will be 
synthesized onboard spacecraft 1. This means that spacecraft 2 and 3 will simply need to 
telemeter to spacecraft 1 the particular combinations of the measurements they have made, 
which enter into the ( combination. Since the space of all the interferometric combinations 
can be constructed by using four generators (a,/3,j, C,) we conclude that spacecraft 2 and 
3 will each have to telemeter to spacecraft 1 four uniquely defined combinations of the 
measurements they have performed. 

The time-delay interferometric combinations require use of phase measurements that are 
time-shifted with enough accuracy to bring the laser phase noise below the secondary noise 
sources. The required time resolution in the time-shifts should be equal to about ~ 50 ns 
for shifts tens of seconds in size. This is because the correct sample of the shifted data 
should be as accurate as the armlength accuracy itself. It can be shown that performing the 
time-shifting, on data sampled at ~ 10 Hz, by using digital interpolation filters, does not 
provide the required accuracy to effectively cancel the laser phase noises (see Appendix B 
for a detailed calculation). 

An alternative approach for achieving a timing accuracy of at least 50 ns would be 
to sample each measurement at ~ 20 MHz or higher and store ~ 2.0 x 10 9 samples in a 
ring buffer for obtaining the data points of this measurement at the needed times. The 
phasemeter would then average these measurements over a fixed time period (perhaps a 
tenth of a second) centered around the sampled times at which the phase measurements are 
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needed. The data is then exchanged among the spacecraft and the TDI combinations are 
formed. 

This method can however be further refined by actually sampling every phase difference 
a few times, each time at ~ 10 Hz, but with a delay between the start time of every sampled 
version of the same phase difference. That is, we envision triggering the phasemeter such 
that the time series are sampled at the times required to form the TDI combinations. In 
this case the limitation of the finite sampling time in the determination of the delayed phase 
measurement is replaced by the timing precision of the phase measurements, which can be 
many orders of magnitude smaller than the smallest sampling time of the phasemeter, as 
we will show below. 

As a concrete example, the («/<,., $<,., 7/ ., Ck>.) basis in the master-slave configuration 
requires measurements {521,821,1,521,3,521,13,531,831,1,831,2,831,12} from spacecraft 1, mea- 
surements {S12, Si2,2, si 2j 3, 812,23} from spacecraft 2, and measurements {si 3 , si 3j2 , 813,3, £13,23} 
from spacecraft 3. By sampling the data s 2 i at the times nj f s , nj f s — L\, n/ f s — L 3 , nj f s — 
L\ — L 3 , where f s ~ 10 Hz is the sampling frequency, and n — 0, 1, 2, ... (and similarly for 
812,813 and S31) we can obtain the entire data at the required times, these being limited 
only by the timing precision of the phasemeters, the time synchronization accuracies of the 
clocks, and the armlength accuracies. This scheme requires sampling each signal four times 
at a sampling frequency (10 Hz) much smaller than what would be needed if sampling the 
data to the granularity required by the TDI combinations. Of course to correctly sample at 
10 Hz we must first ensure that the signal frequency bandwidth is less than 5 Hz to avoid 
aliasing problems. 

In practice, the times at which every sample from a given signal are taken could be 
adjusted every l/f s , in order to protect the quality of the laser noise cancellation against 
drifting armlengths. This requires an adequate model of the spacecraft orbits, which could 
be updated as needed from spacecraft ranging data. 

D. Sampling time jitter 

The sampling times of all the measurements needed for synthesizing the TDI combinations 
will not be constant, due to the intrinsic timing jitters of the digitizing systems (USOs and 
phasemeters). Within the digitizing system, the USO is expected to be the dominant source 
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of time jittering in the sampled data. Presently existing, space qualified, USO can achieve 
an Allan standard deviation of about 1CT 13 for integration times from 1 to 10000 seconds. 
This timing stability translates into a time jitter of about 10~ 13 seconds over a period of 1 
second. A perturbative analysis including the three sampling time jitters due to the three 
clocks shows that any laser phase fluctuations remaining in the four TDI generators will also 
be proportional to the sampling time jitters. Since the latter are approximately four orders of 
magnitude smaller than the armlength and clocks synchronization accuracies derived earlier, 
we conclude that the magnitude of laser noise residual into the TDI combinations due to 
the sampling time jitters can be made negligible. 

E. Data digitization and bit-accuracy requirement 

As shown in figure 5, the maximum of the ratio of the laser noise and of the secondary 
noises phase fluctuation amplitudes occurs at the lower end of the LISA bandwidth, and is 
~ 10 10 at 0.1 mHz. This corresponds to the minimum dynamic range for the phasemeters to 
correctly measure the laser fluctuations and the weaker signals simultaneously. An additional 
safety factor of ~ 10 should be sufficient to avoid saturation if the noises are well described 
by Gaussian statistics. 

In terms of requirements on the digital signal processing subsystem, this dynamic range 
implies that approximately 36 bits are needed when combining the signals in TDI, only to 
bridge the gap between laser frequency noise and the other noises and gravitational wave 
signals. More bits might be necessary to provide enough information to efficiently filter the 
data when extracting weak gravitational wave signals embedded into noise. 

VI. SUMMARY AND CONCLUSIONS 

A comparative analysis of different schemes for implementing Time-Delay Interferometry 
with LISA has been presented. In particular, we have shown that the master-slave configu- 
ration is capable of generating the entire space of interferometric combinations identical to 
that derived by using the one-way scheme. This was done under the assumption that the 
noise from the optical transponders was negligible. Our analysis of the phase-locking control 
systems forming the optical transponders shows that this is a valid assumption, indicating 
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FIG. 5: Phase fluctuations spectra are plotted versus Fourier frequency for: (upper curve) raw laser 
noise having spectral density 2.3 x 10 _1 /~ 8 / 3 + 1.4 x 10~ 9 f~ 27 ^ 5 cycles 2 Hz -1 , and (lower curves) 
residual noises entering into the various TDI combinations. The armlength has been assumed to 
be equal to L = 16.67 sec. 

that the noise introduced can be expected to be ~ 1 /icycle/vldz- 

A comparison of the hardware required for each scheme shows that the subsystems needed 
are almost identical, with the only difference being the number of frequency stabilization 
systems. This difference is perhaps not significant when redundancy options are considered. 
The main disadvantage of the one-way method is that the laser frequencies might be offset 
by several hundred megahertz, given the currently envisioned optical-cavity-based frequency 
stabilization systems. This places challenging constraints on the photoreceiver bandwidth 
and bandwidth stability. On the other hand, the master-slave configuration has no such 
problem, allowing the beat-note on the photoreceiver to be the minimum determined by 
the Doppler shift. However, there may be concerns with the non-local nature of the phase 
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locking system, since its performance could be influenced, for example, by pointing stability 
(which also has implications on lock acquisition). Further studies on these issues should be 
performed. 

Given the similarities between the two schemes, in principle either operational mode 
could be implemented without major implications on the hardware configuration. Ultimately 
detailed engineering studies will identify the preferred approach. 

A derivation of the armlengths and clocks synchronization accuracies, as well as a deter- 
mination of the precision requirement on the sampling time jitter, have also been derived. 
We found that an armlength accuracy of about 16 meters, a synchronization accuracy of 
about 50 ns, and the time jitter due to a presently existing Ultra Stable Oscillator will 
allow the suppression of the frequency fluctuations of the lasers below the level identified 
by the secondary noise sources. A new procedure for sampling the data in such a way to 
avoid the problem of having time shifts that are not integer multiples of the sampling time 
was also presented, addressing one of the concerns about the implementation of Time-Delay 
Interferometry. 
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APPENDIX A: FUNDAMENTAL LIMIT OF THE LISA PHASE TRANSPON- 
DER 

The following is a derivation of the noise added by the optical phase measurement. This 
calculation only considers noise added by the detection system, it does not include the 
20 /^cycles/ vTlz due to the optical-path noise, and it is purely quantum mechanical [19]. 
The reason for performing the calculation in quantum mechanical terms is to highlight that 
the measurement process itself does not add shot noise and that, in principle, a perfect phase 
measurement of the field could be made. 
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Let us consider the optical configuration shown in figure Mjo), which is equivalent to 
the optical arrangement for the LISA interferometer (figure Efa)). Let the annihilation 
operators for the local and distant lasers be a and b respectively. We can represent these 
operators as the sum of an average (complex number) component and an operator component 
representing the field fluctuations of zero mean values: 

a = a + 5a , (Al) 
b = (3 + 5b . (A2) 

In Eqs. (jAll IA2|) we define a and (3 to be real numbers. We will assume throughout these 
calculations that 5d,5b <C a,f3 and so terms that are of second order in these quantities 
will be ignored. This approximation holds even for the low intensities found in the LISA 
interferometer. 
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FIG. 6: (a) Simplified optical arrangement for phase locking in LISA, (b) Equivalent optical layout 
for modeling purposes, e is the intensity reflectivity of the beam splitter. 

In an offset phase locking system the field b is offset from field o by a radial frequency 
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. In this case equations IA"T1 and lA"2l become 



a = a + 5a 

b = f3e iuJbt + 5be iuJbt . 



(A3) 
(A4) 



The annihilation operator for the field at the photoreceiver, c, will contain some fraction of 
a and b, 



c = \/l — e a + y/e b . 



(A5) 



The photoreceiver measures a quantity proportional to the photon number of this field, 



n r = c'c, 



&c 



(1 - e) [a 2 + a(5a + 5a))] + e (3 2 + f3(5b + 5ft) 

+yfeVl zr ~e \a(3(e iuJbt + e~ iulbt ) + a(5be iuJbt + 5fte~ iuJbt ) + /3(5ae~^ 6 * + 5a)e iuJbt ) 
(1 - e) [a 2 + a(5a + 5a))] + e \f3 2 + (3(5b + 5ft) 

(a(3 + a(5b + 5b ] ) + [3(5a + 5a 1 ")) cos(uj b t) 



(A6) 



+ (a(io"6 — i5ft) — (3(i5a — iSa))) sinf^t) 



(A7) 



To simplify the notation, we define the quadrature operators which represent the fluctuations 
in the amplitude (5X + ) and phase (8X~) quadratures of the operators a and b, 



5X- b 
5X~ h 



(5a + 5a^) , 
(5b + 6ft) , 
i(5a — 5a)) , 
i(5b - 5ft) . 



(A8) 
(A9) 
(A10) 
(All) 



Each of these quantities is an observable of unit variance for a coherent state (idealized 
laser). Substituting these expressions into equation IA7I we obtain, 



[l-e] 



a 2 + a5X" n 



+ e 



(3 2 + (35 X b 



+2^VT^1 [(af3 + a5X+ + (35X+) cos(u b t) + (a5X b ~ - p8X~) sm(u b t)] (A12) 

This equation contains three terms which can now be identified. The first two terms arise 
from the intensities of fields a and b respectively. These terms are non-interferometric in 
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nature and contain the intensity fluctuations of the individual input beams scaled by the 
efficiency of coupling to the photoreceiver (beam splitter ratio). The third term represents 
the interference between the two fields and provides a beat note at frequency Ub, the differ- 
ence frequency of the two fields. This beat note itself has two parts, an intensity noise part 
oscillating as cos(c^t), and a phase difference part oscillating as sin^t). 

The phase difference can be obtained, for example, by using a mixer to demodulate the 
beat note down to zero. Mathematically, this is a multiplication by sin^i). Terms with 
a cosine multiplier will only exhibit higher harmonics whereas terms with a sine multiplier 
will mixed down to base band frequencies. 

After low-pass filtering the mixer output, we obtain the following expression for the error 
signal 

V e = Q-^-aSX*^ + \PSX^ + V~eVT^~e(a5X- - (35X~) , (A13) 

where we have ignored the noise of the oscillator used in the mixing process. The quantities 
X+ Ub and are the amplitude quadrature operators evaluated at the offset frequency. 

The factor of 1/2 in the coefficients of these terms arises as we have only taken the sine 
component of the intensity noise at the heterodyne frequency. For heterodyne frequencies of 
10 MHz and greater the intensity noise is shot noise limited thus ((A+^J 2 ) m ((X^ ) 2 ) 1. 
The important point is that the error signal is proportional to the intensity noise of each 
beam and the relative phase noise of the two lasers. 

Assuming perfect phase locking, the error signal V e is driven to zero by actuating on 
the phase of a. Setting V e = we find that the controller will attempt to force the phase 
quadrature fluctuations to be, 



Ultimately, what is of interest is the difference between the phases of the fields d and a. 
The annihilation operator for the outgoing beam, d, is also made up of a linear combination 
of a and b 

d = i/ea -y/l-eb , (A15) 
which implies the following phase quadrature fluctuations of d, 5X^ 

5X~ = Ve5X~ - VT^l5X b - . (A16) 
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Under the operational configuration of perfect phase locking, we can substitute into the 
equation above the expression for 5X~ given in equation IA14I 



S Xi = - ^l)6X; + 2ff + ^ . (A17) 

Since the power of the laser A (proportional to a 2 ) will be a factor of ~ 10 8 larger than 
the power of laser B ( proportional to /3 2 ), and also that {(5X^) 2 ) ~ 1, i.e. the signal laser 
intensity is approximately shot noise limited at the heterodyne frequency, we find that Eq. 
flA17|) can be approximated as follows 



V€dX b + 



(A18) 



The quantity of interest is the phase fluctuations in radians. To compare the phase 
fluctuations between the two fields we need to normalize the phase quadrature operator by 
the square root of the average photon number. For d, this means dividing by \fea 

5^=^, (A19) 

v-lh+W^h- (A20) 

Thus the phase of the incoming beam will differ from the phase of the outgoing beam by 
the following amount 



r « n y/l-e 5X+ 
64> d -6<i> b = p , 

while the root-mean-squared value of the phase error can be written as 



(A21) 



e(l^+J 2 > 



1 4e rib 

Here fib = /3 2 is the average number of photons in the weak signal laser, and (\SX£ \ 2 ) is 
the variance of the local oscillator intensity fluctuations relative to the variance of quantum 
noise. Thus the error depends on two parameters, the beam splitter ratio and the intensity 
noise of the local oscillator. Assuming the local laser intensity is shot noise limited at the 
modulation frequency and a beam splitter ratio of approximately 100:1 (e = 0.99) gives a 
phase error with a standard deviation, cu, of approximately 5% of the shot noise limit or 
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less than 1 /xcycle/\/Hz error. This error is therefore much less than the 20 //cycles/vHz 
optical path noise. 

If necessary, this source of error could be removed by altering the detection system to add 
a second detector allowing subtraction and consequent cancellation of the intensity noise of 
the local oscillator laser. Cancellation factors of 100 are readily achievable, albeit with a 
slight increase in system complexity, effectively removing this error contribution entirely. 



APPENDIX B: INTERPOLATION ERROR FOR GENERATING SHIFTED 
DATA POINTS 

Let g(t) be the true signal. It is sampled at intervals At = l/f s , where f s is the sampling 
frequency, to produce the discrete data g[n] = g(t + n/ f s ), n = —1, 0, 1, .... Consider the 
problem of estimating g(t) at some time which falls in between two sampling times, i.e. at 
time t where < \t — no/ f s \ < At, for hq the value of t/ At rounded to the nearest integer. It 
is well known that for an infinitely long dataset, this estimation can be done without error 

n 

using the Shannon formula (20j, assuming that the signal has zero power above the Nyquist 
frequency (f s /2). The error in the estimation with a finite digital filter can be approximated 
using a truncated version of the Shannon formula. This digital filter might not be the best 
one for all signals, but it should be sufficiently close to the optimal filter. The estimated 
function is given by 

N 

9N{t) = ^2 g[n + no] smc(f s t - n - n ) , (Bl) 

n=-N 

for a digital filter of length 2N + 1, and where sinc(x) = sin(7nr)/7nr. The estimation error 
is e N (t) = g{t)- g N {t). 

Changing the sampling frequency will not improve the function estimator for a fixed 
digital filter size: a lower sampling frequency would be insufficient to represent the high 
frequency signal components, and a higher sampling frequency would reduce the size of the 
interval over which the function is sampled to build the estimator. 



Assuming g[t) to be a wide sense stationary stochastic process 2l( with autocorrelation 



function R g {\t\ —^1) = -R 9 (^i,^2) = E[g{ti)g{t2)} (E\\ denotes the expectation value), it 
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follows that the autocorrelation function of the estimation error is equal to 



Re N (tl,t2) — Rg(ti — t 2 ) + 



N 



) 



n + n\ — m — n 2 



sinc(/ s ti — n — rii) sinc(f s t 2 — m — n 2 ) 



m,n=—N 



A' 



^ R g {ti -n/f s - n 2 /f s ) sinc(/ s t 



2 — n — n 2 ) 



n=~N 



N 



^ # g (*2 -n/f s - n 2 /f s ) sinc(/ s t 



i — n — rii) 



(B2) 



where rii (n 2 ) is the value of t±/At (t 2 /At) rounded to the nearest integer. This equation 
shows that e^it) is not wide sense stationary; however, it is wide sense cyclo-stationary, 
since Re N (t\,t 2 ) = R eN (t 1 +mAt,t 2 +mAt) for every integer m. This is just a consequence 
of the error varying quasi-periodically with the interpolation time; it is zero when t matches 
a sample time, maximum at the middle between two sample times, etc. 

A fair estimate of the estimation error magnitude can be obtained by considering the 
stochastic process eAr(t) = ejv(t + 9), where 6 is a random variable uniformly distributed in 
[0, At]. eAr(t) is wide sense stationary, and its autocorrelation function is [21] 



The Fourier transform of Re N (r) gives an estimate of the spectrum of the noise induced 
by the digital filters interpolation errors. In particular, R Sn (0) is the broadband standard 
deviation of the noise. 

Taking g(t) to be a laser phase noise with a power spectral density that scales like l// 2 , 
and restricting attention to the frequency range 0.1 mHz < / < 1 Hz, one can calculate 
numerically that R BN (0)/R g (0) = 9 x 1(T 7 for N = 10. Therefore, a filter with iV = 10 is 
good enough only to produce a broadband error on the shifted time series that is ~ 3 orders 
of magnitude smaller than the laser phase noise amplitude. Changing TV in the numerical 
integrations shows that Rg N (0) scales roughly like 1/N. This implies that it would be 
impossible to use digital filters on a slowly sampled time series to achieve the levels of noise 




(B3) 
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cancellation required by the TDI combinations. 
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